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Abstract 

We consider the Stochastic Partial Differential Equation 

Ut = kAu + x)u, t > 0, x 6 Z d . 

The potential is assumed to be Gaussian white noise in time, stationary in space. 
We obtain the asymptotics of the almost sure Lyapunov exponent j(k) for the 
solution as k — ► 0. Namely 7(k) ~ i n (i%) > w here the constant Co is determined by 
the correlation function of the potential. 

1 Introduction 

We study the long time asymptotic behavior for the solution of the parabolic equation 
with random potential 

u t = kAu + £,(t,x)u (1) 

t > 0, xez d 
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with the initial data u(0,x) concentrated at the point x = 0. The potential £(t,x) is 
a mean zero stationary Gaussian field, which is white noise in time. The almost sure 
Lyapunov exponent of the solution is defined as the following limit 

. . , In lift, 0) . , 

7 (« = lim )±J. , (2) 

provided the limit exists a.s. and is nonrandom. The existence of the Lyapunov exponent 
for the above type of potentials easily follows from the Kingman's subadditive ergodic 
theorem [|J. Indeed, the fundamental solution q(s,x,t,y) satisfies 



when s < t < u, while 



q(s, 0, u, 0) > q{s, 0, t, 0)q(t, 0, u, 0) 



E In q(t, 0,0,0) In Eg(t, 0,0,0) 



t ~ t 

which is bounded uniformly in t [§]. Thus, due to time stationarity of the potential we 
can apply the subadditive ergodic theorem to the process \nu(t, 0), and thus prove the 
existence of the limit in (|j). We are interested in the asymptotic behavior of the Lyapunov 
exponent j(k) for small k. In |2| the bounds on 7(k) were obtained with the lower bound 
being of order t , j 1 , , , and the upper bound of order ^t^t^) • In @ the upper bound was 
improved to have the same order as the lower bound. In the present paper we obtain the 
asymptotics for the Lyapunov exponent. That is, instead of the upper and lower bounds 
with different constants, we show that a constant Cq exists such that 

We shall prove the following theorem: 

Theorem 1.1 Let (W(t,x),t > 0, x G Z d ) be a mean zero Gaussian field with covariance 
E(W(t, x)W(s, y)) = min(s, t)Q{x — y), where Q(x) is not identically constant. Then 
there exists a constant Cq > such that the solution of the stochastic PDE 

du(t, x) = nAu(t, x)dt + u(t, x) o dW(t, x) (3) 

with initial data u(0, x) = 5(x) satisfies almost surely 

\nu(t,0) c 

hm = 7(/c) ~ — as k — > . 

t^oo t n ' lnl/« 

The outline of the proof is the following: Via Feynman-Kac formula we represent the 
solution u(t, 0) as an integral over the set of paths of a continuous time random walk 
starting at and ending at at time t. The contribution from each path is a functional of 
a Gaussian random variable. We study a small subset of the set of all paths which gives 
the main contribution to the solution. In order to do that we need to study the structure 
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of the maxima of the Gaussian field defined over the set of all paths of the random walk 
with a fixed number of jumps. 

The paper is organized as follows. In section |2| we introduce the notations. In section 
|3] we prove the lower bound for the Lyapunov exponent, while in section ^ we prove the 
matching upper bound. 

2 Notations and Preliminary Considerations 

Let P m be the set of paths of a discrete time random walk of length m starting and finishing 
at 0. Let p(t, m) be the probability that at time t the Poisson process with intensity 2dtt 
is equal to m. Let S(t,m) be the set of possible times of the jumps: t = (t 1 , ...,t m ) with 
< t\ < ... < t m < t. For t G S(t,m) and x = (x ,Xi, ...,x m ) G P m we define 

X m (t,x) = W(t u x )-W(0,x ) + W(t 2j x 1 )-W(t u x l ) + ... + W(t,x m )-W(t m ,x m ) . (4) 

By the Feynman-Kac formula the solution of equation |] can be written as 

00 1 r -ml — 

u(t, 0) = £ p{t, m) £ — - / -e*^d tl ...dt m = 

Jsfc - rm t<=S(t,m) 
oo . 

]T K m e- 2dKt / e Xm(t ' x) dt 1 ...dt m . (5) 

m=0 P ~ 

xfc ^ m t&S(t,m) 

The proof of both lower and upper bounds will be a refinement of the corresponding 
considerations in || and [Q. The central idea of the lower estimation in is the fol- 
lowing: with the fixed times of the jumps (ti,...,t m ) one selects a path, which gives a 
significant contribution to the integral in (|5]). The path is constructed as follows. With 
(xq, Xk-i) selected, one takes x^ to be the neighbor of x^-i, which maximizes the 
increment W(tk+i, Xk) — W(tk, Xk). Such a path follows the local maximum of the poten- 
tial. Better estimations can be obtained by similar optimization with respect to pairs, 
triples, etc. of random jumps. In this paper, rather than maximizing the increments 
W(tk+i,Xk) — W(tk,Xk) at each time step, we shall consider the maximum of the filed 
X m (t,x) over all possible paths x and times of the jumps t. We further find the optimal 
number of jumps as a function of t and n. 

The proof of the upper estimate is similar to that of ||. However, now we are able to 
keep track of the constant, so that to make sure that the upper and lower bounds match, 
giving us the exact asymptotics for the Lyapunov Exponent. 

For a nonnegative let mo(a,t) = [at] (the integer part). Later a will be taken to be 
equal to zln - (1/k) where z will vary. It is the paths with the number of jumps of order 
ln" 2 (l//t)t that give the main contribution to the positive solution u(t,x). Let us stress 
that for typical paths the number of jumps mi has the order nt, that is mo >> mi. 
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Let T a be the space of paths of the continuous time random walk starting and finishing 
at x = with no more than m Q (a, t) jumps. Thus to specify an element of T a we need to 
specify the number of jumps m, t G S(t,m), and x G P m . Then X m (t,x) is a Gaussian 
field over T a . Let 

g(a,t) = sup X m (t,x) , and /(a, t) = Eg (at, t) . (6) 

(m,t,x)&T a 



-OO 



Then f(at, t) is a superadditive nonrandom function of t. In order to see that lim t - 
exists we need to show that / does not grow faster than linearly in t. We shall use the 
following entropy estimate [|]]: For a mean zero Gaussian field X(r) over a set T one 
defines the canonical metric d(ri, r 2 ) = (E(X(ti) — X(r 2 )) 2 ) 1 / 2 . The metric entropy N(e) 
is the smallest number of closed <i-balls of radius e needed to cover T. One has the 
following theorem |jj 

Theorem 2.1 There exists a universal constant K such that 

oo 

EsupX(r) < K j ^\nN(e)ale 
T o 

provided that the RHS is finite. 

Let us estimate the entropy function N(e) for the field X m (t,x) defined over T a . Below 
ci, C2, etc. denote constants, which depend only on the field W(t, x), unless it is indicated 
otherwise. Note that the diameter of T a does not exceed 2y/i, since E(X m (t,x)) 2 = t. 
Assume that e < 2y/i. For t = (t%, ...,t m ) and s = (si, s m ), with m and x G P m fixed, 
the distance between (m, t, x) and (m, s, x) in the canonical metric is estimated as follows: 



d((m,t,x), (m,s,x)) = y/E(X m (t, x) - X m (s, x)) 2 < c lt 



k^2\U- Si | • (7) 

\ i=l 

The inequality is the statement of Lemma 2.1 of 0, it can also be seen directly from (|j). 
Consider the finite subset U(t,m) of S(t,m) defined as follows: t = (tx, ■■■,t m ) G U(t,m) 
if < t\ < ti < ... < t m < t, and all are positive integer multiples of By (|^) 

the elements (m,t,x) with t G U(t,m) form an e-net in T a . The number of elements in 

([ c i tm ] m i); 

[/(£, m) is equal to — |r , which, for e < 2yt is estimated as follows 

([^]-l)!m! 

([^]+m-l)! < < cjtm J_ < 

([^] _ i)! m i ~ m! ~ £ 2 m! ~ £ 2 ' 

where the last inequality is due to Stirling's formula. Since there are less than (2d) m 
elements in P m and m < m , the entropy is estimated as follows for e < 2\fi 

m ° r t r t 

m=0 £ £ 
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Thus, 



Therefore 



N(e) < max{l,(^) m °} for all e 



^ln N(e)de < ^frriQ J ylii(-^-)de = V c^mot J yln(— )cfe < c^y/mot . 

The entropy estimate implies that 

f(a, t) < c 6 y/mo~i = c 6X J[at]t < c 7 (a)t . 

Therefore the limit lim^oo -^y^ exists and is finite. It will be denoted by F(a). Due to 
the fact that W(t,x) is a Wiener process in t, and is therefore self-similar, 

lim -EsupX m (t, x) = \fa lim -EsupX m (t, x) , 

i->oo t T a t Ti 

which implies 

F(a) = V^F(l) . (9) 

In order to relate the maximum of a 'typical' realization of a Gaussian field to the 
expectation of the maximum, controlled by the entropy estimate, one uses the Borell's 
inequality [|l|] 

Theorem 2.2 LetX{r) be a centered Gaussian process with sample paths bounded almost 
surely. Then Esup T X(r) < oo, and for all A > 

Prob{| supX(r) - EsupX(r)| > A} < 2e^ A2 /- 2 ? ( 10 ) 

T T 

where a 2 = sup T EX 2 (r). 

We next make an intuitive argument providing us with the asymptotics for the Lyapunov 
exponent. The rigorous derivation will be given in the following two sections. We claim 
that the main contribution to the sum in the RHS of ([5]) comes from the terms with 
m ~ zlnT 2 {\ / n)t , where z does not depend on k. Moreover, since we are taking the 
logarithm of u, only the factors n m and e Xm<yt,x ' i are important. We substitute the latter 
by e F<yZln t 1 /"))*. Taking the logarithm of the product of these two factors, that is ignoring 
the summation and integration in the RHS of we arrive at the expression 

ln(K m e F(zln ~ 2(1/K))i ) . 

Dividing by t we obtain 

fi_ F{lh _ F(l)Jz-z 



ln(l//«) 



The maximum of this function is equal to 4 and is achieved at z — F(l) 2 /A. We shall 



41n(l/re) 

prove that j(k) ~ 4wi/^) ■ Note that F(l) is greater than zero if Q(x) is not identically 
constant. 
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3 Proof of the Lower Bound 

We select Co (which is the same constant as in the statement of theorem |0|) as follows: 
c = F ty . Let So > be given. In order to prove the estimate from below it is sufficient 
to show that there exists k such that 

7(«) > wlk) (1 " 6o) f ° r K ~ Ko ■ 

Let a(n) = ln a ■ Since we know that the Lyapunov exponent exists and is nonrandom 
in order to prove (|ll]) it is sufficient to show that for k < k , for large enough t 

Prob{iln( £ K m e~ 2dKt J2 J e^^dt^d^) > ~ } > 1/2 . (12) 



m<[a( K )t] x&Pm teS(tm) 



Since > ^jr^ for m < [a(/c)£], the LHS of the first inequality in (|I2) can be 

estimated from below by the following expression 

C o , , 1 



2ck + -ln( sup sup / e Xm{t ' x) dt 1 ...dt m ) . 



ln(l/Kj t m<[a(«)f] £ 6 p m 

teS(t,m) 

Select /«i small enough so that 

2dK < 2lnTIM fOT " " Kl ' (13) 
Thus what we want to show is that for sufficiently small k 

Prob{^ln( sup sup / e^^dt^dtm) > r7 ^(2 - |)} > 1/2 (14) 
i m<[a(K)t]xeP, J m[l/K) I 



for large t. 

We first consider the field X m (t,x) to be defined on T^, a space of paths somewhat 
smaller than T a . Namely is the space of paths of the continuous time random walk 
starting and finishing at x = with no more than mo(a,t) = [at] jumps, with the jumps 
separated from each other and from the endpoints of the interval [0, t] by a distance of at 
least 2r. Let 

g r (a, t) = sup X m (t, x) and f r (a, t) = Eg r (a, t) . 



Then f r (a,t) is a superadditive nonrandom function of t and there exists the limit 
F r (a) = Umt-»oo As in © 
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We next demonstrate that F r (l) — > -F(l) as r — > 0. Indeed, let an arbitrary e > be 
given. Select £ such that F(l) - < e/2. Since / r (l,t ) -> /(Mo) for fixed t , we 

can find r such that -F(l) — - < £ for r < ro- Note that - ^ is increasing in £. 

Therefore -F(l) — lim^oo - ^ < £ for r < r . Thus -F(l) — F r (l) < e for r < r , which 
proves our claim. 

Now we estimate from below the supremum of the field X m (t,x) over T^. 



Let 5± = So ^° . Select k 2 such that 



F a ^(l) > F(l) - -± for k < k 2 . (15) 



Then for k < k 2 there exists ix = ^i(^i, K ), such that for t>t\ 



Jain) Si 

Esup X m (t,x) >t(F\a(K)) v 



t(V^)F a(K) (l) " X^^) > yfi(dt(F{l) - 50 . 



Then by Borell's inequality ( |Tol) with a 2 = t and A = \J a{n)5it 

Prob{sup X m (t,x) < v /a(/c)t(F(l) - 2^)} < 26^ < \ , (16) 

T"l v 4 

where the last inequality holds for t sufficiently large. The fact that we can estimate the 
supremum of the field X m (t,x) from below does not immediately allow us to prove (|I4]) . 
The problem is that the expression in the LHS of the first inequality in (|14D involves 
integration in ti, ...,t m , and thus we need to study the fluctuations of the field X m (t,x) 
as t is varied. Consider the field 

Y{t x ,t 2 ,m,x) = X m (ti,x) - X m (t 2 ,x) 

defined on the set U a> p C T a x T a . We say that ((ti,mi,Xi), (t 2 ,m 2 ,x 2 )) G U a> p if and 
only if mi = m 2 < [at], x\ = x^, and \t\ — t 2 \ < (3 for 1 < i < m\. Let 

G(a, (3) = lim sup -E sup Y . 

Then G(a,P) = y/aG(l, a/3). Below we are going to show that G(l, 0) — » as (3 — > 0. 
First, assuming that we have this statement, we prove (|TTD . We take k 3 so small that 



G(l,a(«)) < 5 1 for k < k 3 . (17) 
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Then for t sufficiently large Esup^^ i Y < 5\ya(K)t. Thus by Borell's inequality with 



er 2 < 2t and A = Ja(n)Sit 



Prob{ sup Y > 2Ja(K)t5 1 } < 2e — * — < 



where the last inequality holds for t sufficiently large. Take kq = min{Ki, K2, K3}, so that 
(0), (H), and dHD hold. 

From dT5[ ) and flTE| ) it follows that for k < k,q for t sufficiently large with probability 
of at least 1/2 one has 



sup X > ya(K,)t(F(l) — 25\) and at the same time sup Y < 2ya(/c)i5 1 . 

T 1 . , U a( K ),l 

Therefore for large t on a set of probability of at least 1/2 

-ln( sup ^sup f e Xm ^dt x ...dt m ) > I l n (e*v / °W( F ( 1 )~ 4<5 i)) = yJ^){F{l)-AS x ). 

t m<[a(K.)t] xSPm ~ •* t 
teS(t,m) 

(19) 

Recall that F(l) = 2y^, ^a(n) = j^L, and 61 = S 0y /c^/8. Thus the RHS of is 
greater or equal than ln /jV) ~~ if) - Thus the estimate (|14D holds. It remains to show 
that 67(1,/?) -> as /? -> 0. 

By (|8|) the entropy function Nx(s) for the field X m (t,x) defined over 7\ (that is with 
m = [<]) is estimated as N x (e) < max{l, The entropy function A^y(e) for the field 

Y(ti,t2,m,x) defined on the set Ui t p is estimated as A r y(e) < N] c (e/2). The diameter of 
U\£ in the canonical metric associated with Y does not exceed 2sup c/ij3 (EF 2 ) 1 / 2 < 2y/]3i. 
Thus for (3 < c by the entropy estimate 



r2W0t , f2W(3t A c + , rvl 3 / c I 

EsupF < K / J\nN Y (e)de < K / \ H— ) n de = 2^2cKt \ Wln(— )<fc . 
v x ,p Jo Jo V £ Jo V £ 

Then, 



1 ^_ z-vwc 1 

G(l,/3) = limsup-Esupy < 2^^ / Wln(— )cfe 
^cxd t u 1>f3 Jo V £ 

which tends to as (3 — > 0. This completes the proof of (11). 



4 Proof of the Upper Bound 



Let 5q > be given. In order to prove the estimate from above it is sufficient to show 
that there exists kq such that 

7(«) < TJTT-^ 1 + ^ for ^ — ^0 • (20) 
ln(l/K) 
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As in section |3] we ust the Feynman-Kac formula for the solution. Since we know that 
the Lyapunov exponent exists and is nonrandom in order to prove (p0|) it is sufficient to 
show that for k < kq for large enough t 

Prob{iln(f; K m e~^ £ / e x ^dU...dt m ) < } > 1/2 . (21) 



m=0 T pp 

x ^ m tES(t,m) 



The m summation in the LHS of the first inequality of (21) will be performed over disjoint 
intervals separately. For 7 > let us estimate the probability of the following event 

A au a 2 (l) = { E K m e~ 2dKt J2 J e x ^dt 1 ...dt m >e^} (22) 

Calculating the number of terms in each of the sums and the area of the domain of 
integration, we estimate the LHS of the inequality in ( p2]) from above by 

a 2 t +m 

t{a 2 -a 1 )\n- 2 (l/K 1 )n ait e- 2dKt {2d)^7^( sup — -)exp( sup X m ) (23) 

a 2 i i TTl. dot -i 

m< m< ry-f- — 

- L ln 2 (l/ft) J -Mn^fl/ftjJ 

Since t(a 2 — ai) In (1/k) < (2d) In2 (V«) the logarithm of ( p3|) is estimated from above by 

2a t t m ~ 

" " ln(2d) - oitln(l/«) - 2ckt + In sup (— -) + sup X m (t,x), (24) 



hi 2 (l/K)~" v '"' " ' ~ m <La* j v ro! m<[ °2« 

— L ln 2 (l/<t) J — L h?(l/K) J 

Let Ci, C2, etc. denote constants which may depend only on the dimension d. From 
Stirling's formula it follows that 

. ,t m C2max(l,a 2 )tln(ln 2 (l//s;)) 
in sup ( — -) < sup m(mt — mm + Ci) < §7 — r\ • 

m<[ ° 2 ' J m! m<[ a ?* J ln (!/«) 

— L ln 2 (l/<t) J — L ]n 2 (l/(t) J 



Therefore the quantity in (|24| ) is estimated from above by 

,, ,~ „ x , . s c 3 max(l,a 2 )tln(ln 2 (l/K)) 

sup X m (t, x) - ait ln(l/«) + V ff- . \ V 7 ;; , 

m<[ j*i ] In (1/k) 

Thus the probability of the event A ai)02 (7) is estimated from above by 

^ , r „ /- ~n , / / \ Csmaxfl, a2)tln(ln 2 (l/«;)) , 
Prob{ sup X m {t, x) - ait In 1/k + 1 if- . \ W " > 7*} = 

m<[ ffl ] 111 (V«) 

^ ,~ ~s , / / s o(/c) maxfl, an) „ , , , 

Prob{ sup X m t, x > t 7 + ai In (1/k - J^-—A-L_^)} . 25 

, r 02* 1 ln(l/K) 
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Recall that Esup <r a 2 t i X m (t,x) is a superadditive function of t with the limit 

EsU Pm<[ ff* ]^m(t,g) fi- 

*-*oc * ln(l/«) V ' 



Thus Esup < r g 2 t ,X m (*,x) < jf^lyF(l). Therefore by Borell's inequality (|10D with 
a 2 = t the probability in fl25|) is estimated from above by 



. 1 . , . . . ^fao „, x o{k) max(l, do) s0 . 
2e X p(-- t ( 7 + 01 KI/k) - ^,.(1) _ M| (26) 

Take 

_ C (1 + ^q) _ c (l + f)-n g 

7 ln(l/«) ' 7n ln(l/«) ' 

where e is a positive number to be selected below. We cover the axis [0, oo) by intervals 

[a", a%] = [(n — l)ei, nsi], n>l, 

where S\ is to be specified below. Then for large t 

oo „ oo 

P(£ K m e~ 2dKt J>2 J e x -^dt 1 ...dt m > e*) < £ P(A a n >a n( Jn )), (27) 

since 

oo 

e 7ni < e 7 * for large t . 

71=1 

Each term in the RHS of (|27| ) is estimated by the expression of the form (f26|). Thus, in 
order to demonstrate that ( pT|) holds it is enough to show that 

V9pW ^fV* i Q i v 7 ^" ei/-| \ o(/t)max(l,Q 2 1 
^2exp(--t( 7 + ^Ty^ - i^n 1 ) } } < 2 (28) 

for large t. Recalling the definition of a™, a 2 , and 7™ we see that 

7 n + a i _ V«2 o(«) max(l, a£) _ 



m(l/«) Ln(l/«) v ' ln(l//s:) 
(cq(1 H ) — ne + nei — 2^/neiCo — E\ — o(«) max(l, raei)) . (29) 



1h(1/k) v v 2 

Since Co + 2 — 2y/c z > for all z the expression in ([29]) is estimated from below by 



Cq— — + n£i — E\ — o(k) max(l, n£i)) 



ln(l/K) v 2 
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Take such kq that o(k) < min(^ a , ~) for k < Kq. Take E\ = and e = y. Then 
the last expression is estimated from below by i^Y/^y n ^if i - Therefore the LHS of (28) is 
estimated from above by 

which can be made arbitrarily small by selecting t large enough. This completes the proof 
of the upper bound. 
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